Spin-Selective Transport of Electron in DNA Double Helix 
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The experiment that the high spin selectivity and the length-dependent spin polarization are 
observed in double-stranded DNA [Science 331, 894 (2011)], is elucidated by considering the com- 
bination of the spin-orbit coupling, the environment-induced dephasing, and the helical symmetry. 
We show that the spin polarization in double-stranded DNA is significant even in the case of weak 
spin-orbit coupling, while no spin polarization appears in single-stranded DNA. Furthermore, the 
underlying physical mechanism and the parameters-dependence of the spin polarization are studied. 

PACS numbers: 87.14.gk, 85.75.-d, 87.15.Pc, 72.25.-b 



Molecular spintronics, by combining molecular elec- 
tronics with spintronics to manipulate the transport of 
electron spins in organic molecular systems, is regarded 
as one of the most promising research fields and is now 
attracting extensive interest [l|-|4| , owing to the long spin 
relaxation time and the flexibility of organic materials. 
Unconventional magnetic properties of molecular systems 
reported in organic spin valves and magnetic tunnel junc- 
tions, are attributed to the hybrid states in the organic- 
magnetic interfaces and to single-molecule magnet 
(H- Organic molecules would not be suitable candidates 
for spin-selective transport because of their nonmagnetic 
properties and weak spin-orbit coupling (SOC) [lOj. 

However, very recently, Gohler et al. reported the spin 
selectivity of photoelectron transmission through self- 
assembled monolayers of double-stranded DNA (dsDNA) 
deposited on gold substrate They found that well- 



organized monolayers of the dsDNA act as very efficient 
spin filters with high spin polarization at room temper- 
ature for long dsDNA, irrespective of the polarization of 
the incident light. The spin filtration efficiency increases 
with increasing length of the dsDNA and contrarily no 
spin polarization could be observed for single-stranded 
DNA (ssDNA). These results were further substantiated 
by direct charge transport measurements of single ds- 
DNA connected between two leads [12[ . Although several 
theoretical models were put forward to investigate the 
spin-selective properties of DNA molecule based on sin- 
gle helical chain-induced Rashba SOC 3, 14 1, the mod- 
els neglect the double helix feature of the dsDNA and 
are somewhat inconsistent with the experimental results 
that the ssDNA could not be a spin filter. Until now the 
underlying physical mechanism remains unclear for high 
spin selectivity observed in the dsDNA ll| 16j . 

In this Letter, a model Hamiltonian, including the 
small environment-induced dephasing, the weak SOC, 
and the helical symmetry, is proposed to investigate the 
quantum spin transport through the ssDNA and dsDNA 
connected to nonmagnetic leads. We interpret the ex- 
perimental results that the electrons transmitted through 
the dsDNA exhibit high spin polarization, the spin filtra- 
tion efficiency will be enhanced by increasing the DNA 



length, and no spin polarization appears for the ssDNA. 
The physical mechanism arises from the combination of 
the dephasing, the SOC, and the helical symmetry. No 
spin polarization could be observed if any aforementioned 
factor is absent. In addition, the spin polarization could 
be considerably enhanced by appropriately increasing the 
dephasing strength or by decreasing the helix angle. 

The charge transport through the dsDNA, illustrated 
in Fig. [TJ can be simulated by the Hamiltonian: 



n = Hdna + n 



lead 



7i c + its 



n d . 
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Here "Hdna = £n=iKf= 
Xnc\ n C2n + H.c] is the Hamiltonian of usual two-leg lad- 
der model including spin degree of freedom 17[ , with N 
the DNA length, ct n = (cL*,ct n ^) the creation operator 
of the spinor at the nth site of the jth chain of the ds- 
DNA, £j n the on-site energy, tj n the intrachain hopping 
integral, and X n the interchain hybridization interaction. 

Ulead + Tic = Z)fe,/3( / 3=L,R)[ e ^fe a ^fe a /3fc + tp<*pk( C lnf, + 

C2na)+H.c] describe the left and right nonmagnetic leads 
and the coupling between the leads and the dsDNA, with 
n L = 1 and riR = N. H ao and Tid are, respectively, the 
Hamiltonians of the SOC term and the dephasing term, 
which will be discussed in the following. 

When a charge is moving under an electrostatic poten- 
tial V, an SOC arises H so — 4m ^ c a VU • (a X p), with the 
electron mass m, the speed of light c, the Pauli matrices 
a = (cr x ,a y ,a z ), and the momentum operator p. In the 
dsDNA, the differences of the potential are usually big- 
ger along the radial direction f than that along the helix 
axis (z-axis in Fig.[I| On the other hand, since the 
differences of V are especially large between the interior 
and the exterior of the dsDNA, dV/dr is very large at 
the boundary r — R with R the radius [l{|. Hence, it 
is reasonable to consider the r-component of V only and 
the SOC can be simplified in the cylindrical coordinate 
system H so = -f a- (r x|) with a = i^r^V^r). Con- 
sidering a charge propagating in one helical chain of the 
dsDNA (e.g., the dotted line in Fig. [T]), the momentum 
p = pnl» with Ih the unit vector along the helical chain di- 
rection. Thus H so is reduced to H so = — ^[ctxP||+P||0'_l], 
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FIG. 1: (color online). Schematic view of the dsDNA with 
radius R, pitch h, helix angle 6, and arc length l a - The circles 
represent the nucleobases, where the full ones assemble one 
helical chain and the open ones form the other helical chain. 
The arc length satisfies l a cos8 = RAtp and l a sin8 — Ah, 
with Atp and Ah being the twist angle and the stacking dis- 
tance between neighboring base pairs, respectively. We set 
R = 1 nm, Ah = 0.34 nm, and Atp = -|, which are typical 
values of B-form DNA. Other parameters can be obtained: 
h = 3.4 nm, 6 « 0.5 rad, and l a ~ 0.71 nm. 



where a±_ (tp) = a x sin ip. smQ—tjy cos ip sin 6+a z cos 9 with 
the helix angle and ip the cylindrical coordinate. Since 
the dsDNA consists of two helical chains, the total SOC 

is H so = £?=]>. 



#£r^]with a { l ] = a ± (tp) 



(2) 

and <j\_ = <7±(ip + it). The interchain SOC has been 
neglected because it is very small due to the potential 
symmetry. By using the second quantization |20j, % so 
can be written as: 



J \rjti) 



T 0') 1„ 



H.c. 



(2) 



(T±(nA(p), and a. 



(2) 
n+l 



where t so = -f x 

<7_i_ (nAip + tt). l a and A(p are, respectively, the arc length 
and the twist angle between successive base pairs. 

On the other hand, a charge transmitting through 
the dsDNA will experience inelastic scattering from the 
phonons due to the fluctuation of each nucleobase around 
its equilibrium position and other inelastic collisions with 
the absorbed counterions in the dsDNA due to the neg- 



atively charged sugar-phosphate backbones [21|. Such 
inelastic scattering will give rise to the lose of the phase 
and spin memory of the charge. To simulate the phase- 
breaking process, Biittiker's virtual lead is introduced by 
connecting to each nucleobase 22|, |23[ , with the Hamil- 
tonian of the dephasing term being: 

H d=Yl ( £ inka] nk a ink + t d a] nk c jn + H.c). (3) 

j,n,k 



a jnk — { a ]nkt> a ]nk{) i s tne creation operator of the vir- 
tual lead and td is the coupling between the nucleobase 
and the virtual lead. 

Let us demonstrate analytically that the ssDNA could 
not behave as a spin filter. In continuous real-space spec- 
trum, the Hamiltonian of the ssDNA containing the SOC 

term is written as H ss = ^ - ^[ojj)|| + + V(Z) 

with V(l) the potential energy of the helical chain. 
By taking a unitary transformation with the operator 
U(l) = e {ima/H 2 ) J^^dl ^ ig transformed into 

U' ss = WH SS U = ll - + V(l), which is indepen- 
dent of spin. Therefore, no spin polarization could be 
observed in the ssDNA, regardless of the SOC term, the 
existence of the dephasing, and other model parameters. 
This result can be obtained also by using the discrete 
Hamiltonians of Eqs. ([1]) and @. Similarly, we can also 
verify that any kind of SOC could not give rise to spin 
polarization in the ssDNA. 

According to the Landauer-Biittiker (LB) formula, the 
current in the oth lead (real or virtual) with spin s can 
be written as (25|: I qs = (e 2 /h) J2 m , s > T qs,ms>(V m - V q ), 
where V q is the voltage in the gth lead and T qSjma > = 
Tr[r 9S G'T ms 'G a ] is the transmission coefficient from the 
771th lead with spin s' to the gth lead with spin s. The 
Green function G r = [G a ]t = [EI - H DNA - H so - 

J2 qs and r 9 s = ~ with E thc incident 

electron energy (or the Fermi energy). is the retarded 
self-energy due to the coupling to the qth lead. For the 
real left/right lead, S L/Rs = -ir L/R /2 = ^Pl/r^l/R' 
while for the virtual leads, Ti r qs = —iT/2 = —iTrpdt d , 
with the dephasing parameter T and PL/R/d being the 
density of state of the leads. Since the net currents 
through the virtual leads are zero, their voltages can 
be calculated from the LB formula by applying an ex- 
ternal bias Vb between the left and right leads with 
Vl = Vb and Vr = 0. Finally, the conductance for spin- 
up (Gf ) and spin-down (Gj,) electrons can be obtained 
G s = {e 2 /h) J2 m s > T-R Sjma >V m /Vh, and the spin polariza- 
tion is P s = (G t - G ; )/(G t + G±). 

For the dsDNA, Sj n is set to e\ n = and 62n — 0.3, tj n 
is taken as t\ n — 0.12 and t2 n = —0.1, and A n = —0.3. 
All these parameters are extracted from first-principles 
calculations 2(| 27 1 and the unit is eV. The helix angle 
and the twist angle are set to 8 = 0.5 rad and Aip — -|, 
which are typical values of B-form DNA. The SOC is 
estimated to t so = 0.01, which is an order of magnitude 
smaller than the intrachain hopping integral. In fact, 
all the results are qualitatively same even for smaller t so - 
For the real leads, the parameters Tl = Tr = 1 are fixed. 
For the virtual leads, the dephasing strength is very small 
with T = 0.005, because the monolayers of the dsDNA 
are rigid and the longest DNA in experiment is short [ll| . 
For this value of T, the phase coherence length is longer 
than the dsDNA and the electron transport through the 
dsDNA will keep its phase coherence [22j. However, the 
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finite dephasing is indispensable. The values of all above- 
mentioned parameters will be used throughout the Letter 
except for specific indication in the figure. 

Figure (3Ja) shows the conductances Gf n and the cor- 
responding spin polarization P s . One notices two trans- 
mission bands — the highest occupied molecular orbital 
(HOMO) and the lowest unoccupied molecular orbital 
(LUMO) — in the energy spectrum, where several trans- 
mission peaks are found for both spin-up and spin-down 
electrons (holes) due to the coherence of the system. For 
the HOMO band, Gf and G4. are almost identical, while 
for the LUMO band, Gf and Gf are very different. A 
bell-shaped configuration is observed in the curve of P s 
vs energy E, where the spin polarization of the dsDNA 
at N = 30 can achieve 0.45, which is comparable with 
the experimental results 



To explore the physical scenario to high spin polariza- 
tion observed in the dsDNA, Figs.EJb) and[2jc) plot the 
conductance Gf and P s in the absence of the dephasing 
(r = 0) and of the helical symmetry (8 = ^), respec- 
tively. It clearly appears that the spin polarization van- 
ishes when F = or 8 = -|, although the conductance is 
still very robust in both cases. When T = 0, the dsDNA 
decouples with the virtual leads and the charge transport 
through the dsDNA is completely coherent. In this case, 
the SOC can not generate any spin polarization due to 
the time-reversal symmetry, the helical symmetry, and 
the phase-locking effect [28[ . A small dephasing is neces- 
sary for the existence of finite spin polarization. Indeed, 
it is reasonable to assume a small T because the dephas- 
ing occurs inevitably in the experiment. On the other 
hand, the spin polarization strongly depends on the DNA 
helix. If there is no helix (8 = ?)j no spin polarization 
could be observed (P s = 0). If the right-handed helical 
dsDNA is transformed into the left-handed one (Z-form 
DNA) with tt - 0, P s (vr -8) = -P s (8) exactly. 

We then focus on the spin polarization P s and the aver- 
aged one (P„), where (P„) = ((G t ) - <G;))/((G t ) + (G+)) 
with (G s ) averaged over the LUMO band. Figures [3Ja) 
and[3]Jb) show P s at a fixed energy E and (P s ) vs length 
N, respectively, for several values of the dephasing pa- 
rameter. One notices that P s and (P s ) are enhanced by 
increasing N at first and then saturate or slightly de- 
cline after a critical length N c . With increasing T, N c 
shrinks monotonically [Fig. EJd)]. The behavior of N c 
vs r can be fitted well by a simple function N c oc r _1 . 
For relatively large T (diamond and triangle symbols), 
P s and (P s ) increase faster in the beginning and satu- 
rate at shorter length with smaller values because the 
device is more open. While for smaller T, P s and (P s ) 
increase slower with increasing N in a wider range of TV 
and have larger saturation values. Let's see T = 0.0004 
for instance (circle symbols). P s and (P s ) will keep ris- 
ing even for N > 100, and P s = 0.34 at N = 40 and 
P s — 0.5 at N = 80. These results are quantitatively con- 
sistent with the experiment [ll[. In fact, the dephasing 
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FIG. 2: (color online), (a) Energy-dependence of conduc- 
tance Gf (solid line), Gj, (dotted line), and spin polarization 
P a (dashed line) for realistic situation, (b) and (c) show Gf 
and P s in the absence of the dephasing and of the helical 
symmetry, respectively. Here N=30. 



has two effects: (i) it promotes the openness of the two- 
terminal device and produces the spin polarization [28| ; 
(ii) it makes the charge lose its phase and spin memories 
and then P s is decreased by further increasing T. Accord- 
ingly, for large T with the phase coherence length [22| 
shorter than the dsDNA length, P s will be quite small. 
P s < 0.05 for T = 0.5 and P s -> if T -> oo. Due to the 
interplay between the above two effects, a small T, rang- 
ing from 0.0002 to 0.01, where is much larger than 
100, is optimal for large P s . In addition, Fig. [3jc) shows 
the averaged conductance (Gf) vs N. (Gf) is declined 
by increasing TV or T, because large N or T will enhance 
the scattering. However, (Gf) remains quite large for 
N = 100 and T = 0.012. Therefore, the dsDNA is a well 
spin filter due to the large P s and (Gf). 

Let us further study the spin polarization by varying 
other model parameters. Figures Ufa) and IUh) show P s 
at E = 0.488 with N = 80 and (P a ) with N = 30, re- 
spectively, as functions of the SOC t so and the dephasing 
strength T. P s and (P s ) are zero exactly when T = or 
i so = 0. Of course, t so is a key factor for the spin polar- 
ization or equivalently < so is "the driving force" of P s . If 
there is no SOC, no spin polarization would appear for 
whatever other parameters are. In general, strong SOC 
usually lead to large (P s ) [Fig. 4(b)]. However, for a fixed 
energy, P s will not increase monotonically with t ao , as 
seen in Fig.QJa). A large P s can be obtained for long ds- 
DNA even for quite small t ao , because the spin polarized 
electrons will accumulate gradually when electrons are 
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FIG. 3: (color online). Length-dependence (a) of P s at E = 
0.488, (b) of (P B ), and (c) of (G t ) for different values of the 
dephasing parameter, (d) The critical length N c vs F. Here 
N c is extracted from the curve of {P a )-N and the solid line is 
the fitting curve with N c oc The legends in (d) are for 

panels (a), (b), and (c). 
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FIG. 4: (color online), (a) P s vs the SOC t so and the dephas- 
ing r at E = 0.488 with N = 80. (b), (c), and (d) show (P s ) 
with A^ = 30 as functions of t BO and T, of T and -, and of ti 
and ti, respectively. 



transmitting along the dsDNA. In addition, we observe a 
large area with red color in Fig.[4]Jb), where (P s ) exceeds 
0.1 for short dsDNA. This implies that the dsDNA would 
be an efficient spin filter in a wide parameters range. 

Figure 0Jc) plots the averaged spin polarization (P s ) 
vs r and ^ by fixing the radius R and the arc length 
l a to account for the rigid sugar-phosphate backbones. 
The helix angle can be changed by stretching the DNA 



molecule [29|. It is obvious that (P s ) is zero in the absence 
of the helical symmetry (8 = -|) and (P s ) is increased 
monotonically by decreasing 8. This indicates that the 
helix of the dsDNA plays a vital role to the existence of 
the spin polarization. Finally, we present the influence of 
the hopping integrals t\ and ti on the spin polarization, 
as illustrated in Fig. HJd). We can see that (P s ) is small 
when t\ and ti have identical sign and become large when 
t\ and t% have opposite sign. Since the sign of the hopping 
integral is sensitive to the type of neighboring nucleobases 
27 1 and to the twist angle Aip [13] , the spin polarization 
could be improved by synthesizing specific DNA molecule 
and putting force along the helix axis of the dsDNA. 

In summary, we propose a model Hamiltonian to sim- 
ulate the quantum spin transport through the dsDNA. 
This two-terminal dsDNA-based device would exhibit 
high spin polarization by considering the SOC, the de- 
phasing, and the helical symmetry, although no spin po- 
larization exists in the ssDNA. The spin polarization in- 
creases with increasing the DNA length. Additionally, 
the spin polarization could be improved by properly mod- 
ifying the hopping integral and decreasing the helix angle. 

This work was supported by China-973 program and 
NSF-China under Grants No. 10974236 and 10821403. 



[1] 



[2] 
[3] 

[4] 



[6] 

[7] 
[8] 
[9] 

[10] 
[11] 

[12] 



Electronic address: sunqf@iphy.ac.cn 

A. R. Rocha, V. M. Gara'a-suarez, S. W. Bailey, C. J. 

Lambert, J. Ferrer, and S. Sanvito, Nature Mater. 4, 335 

(2005). 

A. Fert, Rev. Mod. Phys. 80, 1517 (2008). 

V. A. Dediu, L. E. Hueso, I. Bergenti, and C. Taliani, 
Nature Mater. 8, 707 (2009). 

M. Urdampilleta, S. Klyatskaya, J-P. Cleuziou, M. 
Ruben, and W. Wernsdorfer, Nature Mater. 10, 502 
(2011). 

Z. H. Xiong, D. Wu, Z. V. Vardeny, and J. Shi, Nature 
(London) 427, 821 (2004). 

C. Barraud, P. Seneor, R. Mattana, S. Fusil, K. Bouze- 
houane, C. Deranlot, P. Graziosi, L. Hueso, I. Bergenti, 
V. Dediu, F. Petroff, and A. Fert, Nature Phys. 6, 615 
(2010). 

J. Brede, N. Atodiresei, S. Kuck, P. Lazic, V. Caciuc, Y. 
Morikawa, G. Hoffmann, S. Bliigel, and R. Wiesendan- 
ger, Phys. Rev. Lett. 105, 047204 (2010). 
N. Atodiresei, J. Brede, P. Lazic, V. Caciuc, G. Hoff- 
mann, R. Wiesendanger, and S. Bliigel, Phys. Rev. Lett. 
105, 066601 (2010). 

S. Schmaus, A. Bagrets, Y. Nahas, T. K. Yamada, A. 
Bork, M. Bowen, E. Beaurepaire, F. Evers, and W. 
Wulfhekel, Nature Nanotech. 6, 185 (2011). 

F. Kuemmeth, S. Hani, D. C. Ralph, and P. L. McEuen, 
Nature (London) 452, 448 (2008). 

B. Gohler, V. Hamelbeck, T. Z. Markus, M. Kettner, 

G. F. Hanne, Z. Vager, R. Naaman, and H. Zacharias, 
Science 331, 894 (2011). 

Z. Xie, T. Z. Markus, S. R. Cohen, Z. Vager, R. Gutier- 
rez, and R. Naaman, Nano Lett. 11, 4652 (2011). 



■5 



[13] S. Yeganeh, M. A. Ratner, E. Medina, and V. Mujica, J. 

Chem. Phys. 131, 014707 (2009). 
[14] R. Gutierrez, E. Diaz, R. Naaman, and G. Cuniberti, 

arXiv: 1110.0354 (2011). 
[15] G. L. J. A. Rikken, Science 331, 864 (2011). 
[16] M. Di Ventra and Y. V. Pershin, Nature Nanotech. 6, 

198 (2011). 

[17] G. Cuniberti, E. Macia, A. Rodriguez, and R. A. R6mer, 
in Charge Migration in DNA: Perspectives from Physics, 
Chemistry and Biology, edited by T. Chakraborty 
(Springer- Verlag, Berlin, 2007). 

[18] D. Hochberg, G. Edwards, and T. W. Kephart, Phys. 
Rev. E 55, 3765 (1997). 

[19] In other words, here we consider the SOC induced by the 
boundary-confining potential. All the results are qualita- 
tively same if other kinds of SOC are considered. 

[20] Q.-F. Sun, X. C. Xie, and J. Wang, Phys. Rev. B 77, 
035327 (2008). 

[21] X.-Q. Li and Y.-J. Yan, Appl. Phys. Lett. 79, 2190 
(2001). 

[22] Y. Xing, Q.-F. Sun, and J. Wang, Phys. Rev. B 77, 



115346 (2008). 

[23] H. Jiang, S. Cheng, Q.-F. Sun, and X. C. Xie, Phys. Rev. 

Lett. 103, 036803 (2009). 
[24] Q.-F. Sun, J. Wang, and H. Guo, Phys. Rev. B 71, 

165310 (2005). 

[25] Electronic Transport in Mesoscopic Systems, edited by 
S. Datta (Cambridge University Press, Cambridge, U.K., 
1995). 

[26] Y. J. Yan and H. Y. Zhang, J. Theor. Comp. Chem. 1, 
225 (2002). 

[27] K. Senthilkumar, F. C. Grozema, C. F. Guerra, F. M. 
Bickelhaupt, F. D. Lewis, Y. A. Berlin, M. A. Ratner, 
and L. D. A. Siebbeles, J. Am. Chem. Soc. 127, 14894 

(2005) . 

[28] Q.-F. Sun and X. C. Xie, Phys. Rev. B 71, 155321 (2005). 
[29] J. Gore, Z. Bryant, M. Nollmann, M. U. Le, N. R. Coz- 
zarelli, and C. Bustamante, Nature (London) 442, 836 

(2006) . 

[30] R. G. Endres, D. L. Cox, and R. R. P. Singh, Rev. Mod. 
Phys. 76, 195 (2004). 



